arXiv:1501.03081vl [physics.data-an] 13 Jan 2015 


Math. Model. Nat. Phenom. 
Vol. 10, No. 1, 2015, pp. 10-20 


Noise-Produced Patterns in Images Constructed 
from Magnetic Flux Leakage Data 

Denis S. Goldobin® C Q, Anastasiya V. Pimenova®, Jeremy Levesley 6 , 
Peter Elkington d , Mark Bacciarelli d 

“ Institute of Continuous Media Mechanics, UB RAS, Perm 614013, Russia 
b Department of Mathematics, University of Leicester, Leicester LEI 7RH, UK 
c Department of Theoretical Physics, Perm State University, Perm 614990, Russia 
d Weatherford, East Leake, Loughborough LEI2 6JX, UK 


Abstract. Magnetic flux leakage measurements help identify the position, size and shape of 
corrosion-related defects in steel casings used to protect boreholes drilled into oil and gas reser¬ 
voirs. Images constructed from magnetic flux leakage data contain patterns related to noise inher¬ 
ent in the method. We investigate the patterns and their scaling properties for the case of delta- 
correlated input noise, and consider the implications for the method’s ability to resolve defects. 
The analytical evaluation of the noise-produced patterns is made possible by model reduction facil¬ 
itated by large-scale approximation. With appropriate modification, the approach can be employed 
to analyze noise-produced patterns in other situations where the data of interest are not measured 
directly, but are related to the measured data by a complex linear transform involving integrations 
with respect to spatial coordinates. 
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1. Introduction 

A common challenge in engineering is the need to infer properties of interest from indirect mea¬ 
surements. The direct measurement of the property of interest may be possible but expensive or 
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Figure 1: Ferromagnetic layer of non-uniform thickness in the field H and the coordinate frame 

inefficient compared to an indirect measurement—an example being the determination of methane 
hydrate saturation profiles in seafloor sediments, which can be done efficiently and cheaply by 
indirect means by measuring water salinity profiles from which saturations can be determined 
(e.g., see [QQ |2j). The challenge is particularly widespread in methods employing acoustic and 
electromagnetic wave scattering (e.g., see (3l[4i|5]|), and in the deformation of static fields due to a 
medium’s susceptibility to these fields. An example of the latter is the magnetic flux leakage (MFL) 
method used to assess casing integrity in boreholes drilled into rock formations for the exploration 
and evaluation of oil and gas resources (e.g., see AS H Ell)- In these cases the profiles of interest 
are complex transforms of the measured data involving integrations with respect to spatial coor¬ 
dinates. These transforms can produce non-trivial and unexpected responses in the reconstructed 
(output) data due to noise inherent to the measurement. Even for uncorrelated input data the noise- 
produced patterns in reconstructed data can be correlated, smooth and indistinguishable from the 
actual patterns in an idealistic noise-free situation. (In some cases, the appearance of smooth, 
noise-produced patterns can be treated as a result of noise accumulation in the reconstructed data.) 
Therefore noise-produced patterns cannot simply be filtered out by means of some data process¬ 
ing procedure (which would be the case for uncorrelated noise). Such noise-induced patterns can 
significantly influence the utility of the reconstruction procedure, and knowledge of the properties 
of these patterns is therefore important, and serves as an valuable guide in the development of data 
analysis algorithms and in future device design. 

In this paper we suggest and implement an approach to comprehensively analyse noise-produced 
patterns in MFL-reconstructed data for different device designs. The problem set-up corresponds 
to measurements made with modem devices designed for the MFL inspection of wellbore casings 
described for example in Q. Our approach is based on the mathematical methodology developed 
in [jU within the framework of large-scale approximation, which is relevant for the case of corro¬ 
sive damage. Such model reduction provides the opportunity for a comprehensive analysis whose 
results remain qualitatively valid for the general case beyond large-scale approximation. In the 
analysis of inaccuracies, it is generally sufficient to have estimates of errors. A similar approach 
can be employed for various data reconstruction tasks. 

The paper is organized as follows. In Section 2 we introduce the basic mathematical model we 
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use. In Section 3 we evaluate the patterns in the reconstructed data produced by ^-correlated noise 
in the measured magnetic field data. In Section 4 we analyse the scaling properties of the noise- 
produced patterns for various device designs, compare the cases of different designs, and discuss 
the role of the defect geometry for their resolvability against the background of noise-produced 
patterns. Conclusions are drawn in Section 5. 

2. Magnetic field in the presence of large-scale defects 

The system under consideration is shown in Figure Q} it comprises a ferromagnetic layer subject 
to an external magnetic field parallel to the layer middle-plane; the x-axis is directed along T~L, 
and the z-axis is perpendicular to the layer middle-plane. Within the framework of large-scale 
approximation, the magnetic field within the ferromagnetic layer is nearly homogeneous along 
the x-axis and is parallel to the layer-middle plane, and thus can be represented by a gradient of 
potential (f>(x,y ), which is a function of two coordinates (x,y) only: H ^ = —V(j>(x,y). To the 
leading order, this potential obeys the equation [|8ll 


V 2 -(e-^)V 2 0(x,y))=O, 


( 2 . 1 ) 


where subscript ‘2’ indicates that V 2 operates in the two-dimensional space of (x, y ). The function 
a(x, y) describes the relative loss of the ferromagnetic material; the layer thickness profile 


2C(i,») = 2Coe^ ( ™>, 


where 2( 0 is the thickness of the undamaged layer. 

Outside of the layer the magnetic field is the 3-d gradient of a harmonic potential; 


H {e ^ = -Vrtx,y,z), 


( 2 . 2 ) 


Aif(x,y,z) = 0, 

with the boundary condition on the ferromagnet surface 


(2.3) 


ip(x 1 y 1 z = ((x,y)) = <j>(x,y) 


and (—'V(^) 2 _,.±oo = 'H ■ Within the framework of large-scale approximation, to the leading order, 
one can translate the boundary condition for p from the layer surface to z = 0: 


<p(z,y, 0 ) = (j>{x,y) ■ 


(2.4) 


By virtue of Eq. (12.31) . the potential p of the external field can be represented in the basis of 
exponential functions 



(2.5) 
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Here A(k x , k y )e~ k ^ is the Fourier transform of the deformation of the magnetic potential field, 
ip + Hx, measured on the plane elevated by (z m — ( 0 ) above the layer surface. While the field 
<p(x, y, z m ) can be found from measurements of II 1 ,' Xt \x, y, z m ) or/and H y Xt \x, y , z m ) by plain 
integration with respect to x or y, evaluation of ip(x, y, z m ) from measurements of Il-f )xt \x , y, z m ) 
is non-trivial; the relationship between ip(x, y. z m ) and H < f yA ' ) (x, y, z m ) is non-local. For this rela¬ 
tionship, one can rely on Eq. (12.51) and its ^-derivative; 

(x, y, z) = IJ d k x d k y kA(k x , k y ) e~ k ^ , k = jk* + k\ . (2.6) 


Thus, one can introduce a non-local operator C which transforms ip(x, y, z m ) (a function of two 
variables (x, y)) into (x, y, z m ) (another function of two variables (x, y)): 

d 

- = C(p(x,y,z m ). (2.7) 

This transform is correctly defined for any harmonic function vanishing for z —> ± 00 , and in the 
Fourier space it is a multiplication by k = \k\. For trivial conditions on the boundaries limiting 
the system in the x- and y-dircctions, the inverse operator C ~ 1 is well defined; the operator £ is a 
homeomorphism. 

The aim of our work is to establish smooth large-scale patterns in reconstructed data pro¬ 
duced by the noise in measurements. Hence, the approximate calculation of these patterns will 
suffice, and we may make a further model reduction for Eq. (12.11) and the system under con¬ 
sideration. Specifically, we make this reduction for z m small compared to the scale of defects, 
4>(x, y) ~ <p(x, y, z m ). Hence, Eq. (12.11) yields an equation for calculation of a(x, y) from measur¬ 
able derivatives of <j)(x, y): V20 ■ V 2 CT = A 2 </>. To the leading order, II (ext> (z m ) ~ TL; therefore, 
Eq. (12.11) can be approximately rewritten as 



-A 2 tp(x,y,z m ) 


<9#i ext) 

dz 


( 2 . 8 ) 


3. Delta-correlated noise in measurements 

Let us consider noise in the measurement data; 

f(x,y) = fo(x,y) + a£(x,y), (3.1) 

where f(x, y) is the measured field, fo(x, y) is the actual field without distortion by noise in mea¬ 
surements, and a£(x,y) is the measurement noise with amplitude a and £(x,y) noimalized as 
specified below in the text. We assume that the noise at the neighbouring measurement points 
to be uncorrelated (bias in measurements can be removed by calibration of the measuring sen¬ 
sor), that is £(x,y) can be assumed to be 3-corrclatcd. We choose the following normalization: 
(£ (x,y)£(x\ y’)) =4 6(x — x’) 5{y — y’). Here and hereafter, (...) means averaging over the noise 
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realizations. The amplitude a of normalized noise is determined by the variance of the measure¬ 
ment error £j mCds) = fd neas ) — / 0 ; 

«4'—y> = a ^) = 4 (3. 2) 

x v y 

where h x and h y are the distances between the measurement points in the x- and ^-directions, 
respectively. 


3.1. Measuring f(x, y ) = dHjdx 

Let us first consider the case of measuring f(x, y ) = dH^’/dx, which is routinely performed by 
devices in common usage 0. According to Eq. (12.81) . the reconstructed profile 


cr(x,y) 


^ f ix’C I" d x"f(x",y) 

1 „ r x r x ’ 

—C I d x' J d x"f(x",y). 


(3.3) 


Here one can change the order of operation of the integration with respect to x and the operator 
£, which is differentiation with respect to z. Hence, the purely noise-produced pattern in the 
reconstructed profile a(x , y) is 

£*(x,y) = 


dx' / dx" £f(x",y) 


(3.4) 


In contrast to the noise in f(x, y), which is 5-correlated, the noise-produced pattern £ a (x 1 y) is 
a result of two integrations and the action of non-local operator C ; therefore, this pattern can have a 
smooth component with a finite correlation length which will not be easily recognisable against the 
background of the actual noise-free profile a 0 (x , y). In Figure [2b one can see the profile e a {x, y) 
and its smooth component calculated for f(x,y) = a^(x,y) plotted in Figure [2b. We need to 
evaluate the characteristic magnitude of these smoothed profiles on the scale L x x L y . 

To calculate the characteristic magnitude of e a , one can consider the average value of e a over 
the area S = L x x L y ; 


£a(x,y) 


L X Ly 


£a(x,y) dx dy . 


(3.5) 


As £(x, y) is 5-correlated, its integral is a Gaussian random variable by virtue of the central limit 

theorem. It is more convenient to evaluate the magnitude of £~ 1 e cr (x, y) and then assess how C 
influences the magnitude of its ( L x x L y )-scale component. According to the central limit theorem, 
the area-average 

C~ l £ a (x ) y) = —II jC~ 1 £ a (x,y)dxdy (3.6) 

J-ix^y J JS 
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is a Gaussian random variable with variance 


T2 T 2 0/2 
^x ^y 


dx 1 / dyi / dx\ / dx'/s^x'/.y,) 

0 Jo Jo Jo 

f Lx f Ly r xa r' 2 

x dx 2 dy 2 / dx' 2 / dx 2 £/(x 2 ,y 2 ) 


rxi 


LILIW 


dxi / dyi / dx' x / dx'/ 


4a 2 


T 2 T 2 0/2 

-Gr ™ 


x / dx 2 / dy 2 / dx 2 / dx 2 (£/(x", yO c/(x", y 2 )) 


dxi / dyi / dx' x / dx 7 / 


x / dx 2 / dy 2 / dx 2 / dx" 5(x" - x") %i - y 2 ) 


4a 2 




r»X2 


T2 T 2 0/2 
^ /I 


dxi / dx' x / dx" / dx 2 / dx 2 / dx 2 <5(x , 1 '— x 2 ) 


4a 2 L* 
L 2 x L 2 y U 2 20 


(3.7) 


According to Eq. (12.61) . the action of the operator £ in Fourier space is multiplication by the 
wavenumber k = {k 2 x + k 2 ) 1 / 2 . Hence, for the mode with scale £ x x —i.e., k x = 7r /L x and 
k y = 7T /L y —action of £ is equivalent to multiplication by the factor tt(L~ 2 + Ly 2 ) 1 / 2 . Let a a be 
the characteristic magnitude of the noise-produced profile on the scale L x x L y . Then Eq. (13.71) 
yields the representation 

= \fJJJf) * “ 1 U 


VhJiyfsLy 


L l+ L l 


7r \/(( £ / meas) ) 2 ) h * by 1 L x 
2V5W. 7^4 



(3.8) 


Here we explicitly distinguish the factor 1 /\/~S which features the “reference” convergence law for 
the average value of a noisy variable with respect to the length scale L x , L y . One can observe not 
simply a slower convergence of the amplitude of the noise-produced pattern but its independence 
on the scale length for L x ~ L y . The error is the same on all spatial scales, suggesting that 
measurement ( dH z /dx ) is equally applicable for detection of defects on different spatial scales. 
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Figure 2: Sample realization of noise (a/1-L)£(x, y) (a), and corresponding reconstructed profiles a 
(left column), and profiles a with small-scale oscillations filtered-out (right column), (b): f(x, y ) = 
(dH z /dx ), (c): f(x, y) = H z , (d): f(x, y) — H x , (e): f(x, y) = H y . The smoothed profiles are a 
consequence of noise and should not be misinterpreted as real defects. 
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3.2. Measuring f(x, y ) = H z , H x or H y 


Alternatively to the previous case, some devices measure H z = — dp/dz , H x = —dp/dx and/or 
H y = —dp/dy. In Figured one can see sample profiles £„(x. y) for the cases of measuring 
different components of H. Similarly to Eq. (13.41) for the case of / = dH z /dx, Eq. (12.81) yields 


el Hz) {x,y) = -^C j d x' ef z \x\y ), 

el Hx) (x, y) = f dx ' J dx" e f x) (■ x> \ y), 
£ { a Hy) ( x ,y) = J dx' J d y's { f Hy \x',y'). 


(3.9) 

(3.10) 

(3.11) 


Area-averages of e^ z \ e^a and ' are Gaussian random variables, the variances of which are 
related to the following variances: 



1 




4 a 2 L y L\ 

T 2 T 2 XJ 2 9 

^y ^ 


(3.12) 


r- L 


LI LIU 2 J o 

r% 2 


dx\ / dyi / dx' / dx" 


dx 2 J dy 2 J d x' 2 J dx^ (ef x \x'[, j/i) ef^ix'^ y 2 )J = 


4 ° 2 ^ (3.13) 


LlL 2 y U 2 20 


rL x 


rxi 


ry l 


T2 T 2 0/2 / 
Li x Lty rt J o 

rx 2 


(L/y / dyi / dx' / dyi 


x 


'0 


/ V 2 , . 2 [3 r3 

dj/2 (4^ = L 2 L 2 U 2 V Q X 


Finally, the characteristic magnitudes of patterns on the scale L x x L y are 

V<(4 meaS) ) 2 )^^ 1 L x ^+U y 




,(»») 


y/ZU 

VSLy 

Lx 

K\/WU) K 

Ly l L x L* 

\+ L l 

2V5H 

VsL y 

L l ’ 

^v'«4”'“ l ) 2 > h - 

h v 1 /A* 

+ M 

3H 

VS Uv 

aJ ■ 


(3.14) 


(3.15) 


(3.16) 


(3.17) 
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Table 1: Dependence of the form-factor on geometry of defects 



l{dH z /dx) 

1(H Z ) 

7 (H x ) 

7 {h v ) 

L X ~ Ly 

71 

VTo 

4 

71 

7! 

27r 

Y 

L X Ly 

71 

fM 1/2 

71 

7 r L x 

7T Ly 

2^5 

UJ 

V3 

2V5 Ly 

3 L x 

L X Ly 

7T 

(M 3/2 

7T L x 

-) 

2tt L x 

2a/5 

UJ 

\/3 Ly 

2-y/E Ly 

3 Ly 


It is noteworthy that, in contrast to the case of / = dH z /dx, Eqs. (13. 151) — (13.17b demonstrate the 
decay of the noise-produced patterns on large scale according to the “reference” law oc i /Vs. 


4. Discussion: Scaling of noise-produced patterns, form-factor, 
and optimization of measurement techniques 


In this section we compare 4 possible cases of / = (dH z /dx), H z , H x and H y , and discuss the 
role of the shape of defects for their resolution. 

Let us consider the discrepancy between the scaling laws of amplitude of the noise-produced 
patterns a!'a Hz ^ dl ) and cY ’ 1 : for L x ~ L y ~ L, amplitude (irj >n ^ dx) cx L° is independent of L, 
while (la 1 ' 1 oc L _1 . For the former case, the profile reconstruction technique is equally accurate 
for recognition of both small- and large-scale defects as the noise-produced patterns are of the 
same magnitude on all scales. One can rely on this technique, given the noise in the input signal 
(measurements) can be kept small enough. For the case of / = Hj, the noise-produced patterns 
decay for longer scales and, therefore, the reconstruction procedure is more accurate. Even for 
strong noise the recognition of sufficiently broad defects can be accurate. However, simultaneously 
with good resolution of large-scale defects, reliable identification and characterization of narrow 
defects become problematic or even impossible (it is impossible as well within the framework of 
the virgin equation system, without large-scale approximation). Comparing Eqs. (13.81) and (13.151) . 
one can find the “reference” scale L x *\ measuring / = (dll z /dx) is preferable for defects small 
compared to this scale, while for a larger scale of defects the reconstruction is more accurate with 
f = H z - 


L 


X * 


2\/5 


v / 3\ (( 


((4”S) 2 > 


(meas) \ 2 \ 
b f=dH z /dx) / 


(4.1) 


Accuracy of the reconstruction procedure is influenced not only by defect scale but also by its 


9 



























D.S. Goldobin, A.V. Pimenova, et al. 


MFL: Noise-Produced Patterns 


shape. One can rewrite Eqs. (13.81) and (13.15D — (13.171) as 


with form-factor 7 


a (dH z /d,) = 1{dHz/dx) 


<(e 


(meas; \ 2 


) hx hy 


n 


al Hj) = l( Hj ) 




nVs 


7 {dH z /dx) 


7T L x 
2V5 Ly 


n L :r 


_L y 

Ly L x ’ 




7(^-) - l/ 1 + ^2 > 


7(^) 


\/3 

vr Fz 

2V5 4 

7T 


— + 1 


q \ T L T 

O \ J-Jn 1 J-Jnr. 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 


In Table Q] the behaviour of the form-factor for defects with similar dimensions in the x- and y- 
directions ( L x ~ L y ), defects stretched along the y-axis (L x <C L y ) and along the x-axis ( L x 
L y ) is summarized. Notice, for L y L x , the form-factor is diminished for / = (dll z fdx) and 
/ = and increased for f = H y . For L x L y , the form-factor is increased for all cases of /. 


5. Conclusion 

We have developed an approach for the analytical calculation of noise-induced patterns produced 
when employing a complex non-local transformation for data reconstruction. This approach has 
been applied to the magnetic flux leakage (MFL) method for inspection of wellbore casing in¬ 
tegrity; with the MFL method, the ferromagnetic layer thickness profile is reconstructed from 
measurements of the magnetic field above the layer. Within the framework of large-scale approx¬ 
imation, Eqs. (13.81) . (13.151) — (13.17b have been derived for the device designs based on measuring 
(dH z /dx), II z , H x and/or H y , respectively; these equations form one of the principle results of 
our work. In particular, the case of ( dH z /dx) has been found to be equally effective for all scales 
of defect, while use of H z , II X and II y provide more reliable resolution of broader defects, but are 
inaccurate for the detection of narrow defects. The procedures of reconstruction of the ferromag¬ 
netic layer thickness from H z , H x and II y are not equally sensitive to noise. More importantly, for 
these three cases, the sensitivity depends differently on the defect shape. This dependence is char¬ 
acterised by the form-factor, which, as one can see from Eqs. (14.5b — (14.7b . depends on the long-axis 
orientation for scar-shape defects. Thus, the combined usage of these procedures can be beneficial 
for the accuracy of recognition of defects. 
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